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ABSTRACT 

Context. The planets Uranus and Neptune with small apparent diameters are primary calibration standards. 

Aims. We investigate their variability at 90 GHz using archived data taken at the IRAM 30m telescope during the 20 years 

period 1985 to 2005. 

Metliods. We calibrate the planetary observations against non-variable secondary standards (NGC 7027, NGC 7538, W30H, 
K3-50A) observed almost simultaneously. 

Results. Between 1985 and 2005, the viewing angle of Uranus changed from south-pole to equatorial. We find that the disk 
brightness temperature declines by almost 10% (~ 2a) over this time span indicating that the south-pole region is significantly 
brighter than average. Our finding is consistent with recent long-term radio observations at 8.6 GHz by Klein & Hofstadter (2006). 
, Both data sets do moreover show a rapid decrease of the Uranus brightness temperature during the year 1993, indicating a 

■ temporal, planetary scale change. We do not find indications for a variation of Neptune's brightness temperature at the 8% 
' level. 

, Conclusions. If Uranus is to be used as calibration source, and if accuracies better than 10% are required, the Uranus sub-earth 

■ point latitude needs to be taken into account. 

1 Key words. Planets - Uranus - Neptune 



1. Introduction 

At mm- wavelengths, the planets Uranus and Neptune 
with small apparent diameter are frequently used for cal- 
ibration of astronomical sources and telescope parame- 
ters. In this context it is often tacitly assumed that they 
are constant radiators, at least at short term intervals. 
Observations over a decade or a longer period, preferably 
with the same telescope and receivers and traceable mod- 
ifications, may reveal a long-term variability of Uranus 
and Neptune. 

For a long time, Uranus and Neptune have been sep- 
arated by a distance of ~ 22° or less in t he sky and thus 
allow precise comparative measurements. iGrifBn &: Orton 



( 19931 ) observed Uranus and Neptune in 1990 and 1992 to 
derive their brightness temperature in the millimeter and 



submillimeter regime. However, recent microwave as well 
as observations at visible and near-infrared wavelengths 
indicate that Uranus is variable on tim e scales of several 
years (see e.g. iKlein fc Hofstadteilliooii KH06 in the fol- 
lowing). Due to Uranus' large obliquity of 82°, measure- 
ments from the earth alternate between observations of 
the poles and the equator, i.e. the sub-earth point (SEP) 
latitude varies. KH06 argue that the variations they ob- 
served at 3.5 cm (8.6 GHz) are partly caused by the geo- 
metrical effect, but partly also by temporal variations deep 
inside the Uranus troposphere. However, this variability 
has not yet been studied at millimeter wavelengths. 

We investigate whether these long-term variations de- 
scribed above are noticeable in 20 years pointing obser- 
vations made with the IRAM 30m telescope at 90 GHz. 
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^ This collection of pointing measurements provided also the 
basis for several compilations of quasar flux densities published 
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By relating the observations of the planets to nearly si- 
multaneous observations of the constant secondary cali- 
brators NGC7027, NGC7538, W30H and K3-50A, the 
measurements of the planets are free from changes of the 
telescope performance (reflector adjustments, receiver up- 
grades, etc.) but contain, on the other hand, the errors of 
the secondary measurements. The secondaries are small 
relative to the half power beam width (HPBW) of 27" at 
90 GHz. 



2. Observations, selection criteria 

We use heterodyne observations (0.5 and 1 GHz band- 
width) made between 1985 and 2005 with the IRAM 30m 
telescope (Pico Veleta, Spain) at 86-90GHz (3.4mm). 
Most of the data were obtained during pointing measure- 
ments, or extended pointing sessions during the earlier 
years (see 



Greve et al. 19961). The data arc scans across 



the planets, and other sources, including the fitted Gauss- 
ian profiles, their halfwidths, their pointing offsets, and 
their peak antenna temperature TJ, separately deter- 
mined for the azimuth (Az) and elevation (El) directions. 
An archived measurement is accepted for this analysis if 
the Az and El pointing offsets do not exceed 3-5", i.e. 
being small compared to the HPBW, and if the Az and El 
full width at half maximum (FWHM) are within 3 - 5" 
of the source convolved value. Under these conditions the 
antenna temperatures obtained from the Az and El scans 
agree within ~ 10%, and their average value is used. 

The hot-cold-sky calibration method used at the 30m 
telescope corrects for atmospheric attenuation and gives 
the antenna temperature [K] of the beam convolved 
source. 

The aperture efficiency Cap and fo rward efficiency F^g 
(;see lDowneslll989t iGreve et al.lll998al) are regularly deter- 



mined (i^off to within ± 5 %, Cap to within ± 10 %) and 
the flux density of a point source per beam is Sh = 
2(kM)TXFeff/eap = 3.904 T^Feff/eap [Jy] (with A the 
geometrical area of the 30m reflector, and k the Boltzmann 
constant). It is impossible to recover the actual value 
^off/cap for a certain day in order to derive the flux density 
5b from the archived antenna temperature T^. However, 
from simultaneous observations of the constant sources 
NGC 7027, NGC 7538, W30H and K3-50A (see for in- 
stanc e Steppe et al.lll993 : Reuter fc Kramei]ll998 : Sandelll 
1994I ). with flux densities given in Table 1, we derived the 
gain 5b,soc/JA sec = 3.904 i^cff/cap which we applied to 
the measurements of the planets 5b,pia = pia/^A sec ^ 
5b,soc- The derived flux densities of the planets 5b,pia con- 
tain also the errors of the measurements of the standard 
sources. 

We assume that any time variance of the four sec- 
ondary calibrators is negligible. NGC 7538, K3-50A, and 
W30H are ultra-compact Hii regions which are sta- 



ble. Due t o the flat spectrum of the pla netary nebula 
NGC 7027 (jSanchez Gontreras et"al1 Il998l ) the reported 
flux density varia tion at 1.4 GHz of 0.24% per year 
(jPerlev et al.ll2006l ) IS expected to occur, approximately, 
also at mm- wavelengths. This amounts to a change of 
~ 5 % in 20 years. We note that we usually observed sev- 
eral secondaries per day, using their average to calculate 
the planetary fluxes, thereby reducing this slight effect. 

One pointing measurement consists of 2 Az and 2 El 
scans. The error of a pointing measurement, used in the 
following, is the rms-value of the four scans. Almost all 
observations were made in the elevation range ^ 20° - 
70°. In this range, the gain-el evation correction a t 90 GHz 
is less than 2 % (cf. Table 2 in lGreve et aLlll998bl ). We did 
not correct for this minor effect. 

This procedure does not allow the derivation of abso- 
lute flux densities of the planets since their flux densities 
were used to derive the flux densities of the standards. In 
the following we therefore study relative changes of the 
planet's flux densities and disk- averaged brightness tem- 
peratures. 

We show in Fig.[l] the measurements of the four sec- 
ondary calibrators over the entire observing period and 
list their mean temperatures in Table[TJ An accuracy of 
better than 10% is achieved for all secondaries illustrating 
the precision of repeated pointing measurements made by 
many visiting observers usually for their own programs. 
The observed scatter of antenna temperatures includes 
possible changes of the telescope efficiencies from day to 
day and over the years. Such changes are calibrated out 
when studying relative temperature variations measured 
on the same day. 

From the relation between the observed flux density 
5b, the Planck function B at temperature T{9,tp), the 
beam pattern P{9, and the solid angle U of the planet 
subtended at the time of observation 



5b =/ B[T{9-9\ij~ij')]P{9',ij')dfl' (1) 

"'0 

we obtain when assuming a constant brightness tem- 
perature T{9,ip) ~ Tb H(0, t/j) across the disk (H) of the 
planet 

5b = (2k/A2)TRj / U{9-9',tl^-ij')Pi9',ij')dn' (2) 
Jo 

and with the Rayleigh-Jeans temperature 
Trj = hi./k(exp(hi./(krB)) - i)-i 

-1 



Trj - (AV2k)5b l^j^npdn' 

Trj = {X' /2k) S^{nK)-\ 



(3) 
(4) 



by ISteppe et al] (|l988l . 1 19921 . Il993l ) and iReuter et al] (|l997h . 

Unfortunately, most of the data of 1988 and 1989 are lost, for 
unknown reason. 



The first-order Planck correction at 90 GHz is 
Trj « Tb - hi^/(2k) = Tb - 2.1 K. 

In the above equations, we use for the b eam pattern 
P(M') a Gaussian profile of 27" HPBW jCreve et alj 
1998al) . as it remained constant throughout the years by 
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Fig. 1. Observations of the secondary calibrators at 
90 GHz. Dashed hnes show their rms scatter (cf . Table[T]) . 



using receivers of similar illumination taper. Using the 
IRAM in-house program planets we derive from the flux 
densities the corresponding disk-averaged brightness tem- 
perature Trj [K]. This program takes into account the 
apparent diameter of the planets 6's and the HBPW. In 
Eq.m the correction fact or K for a non-pointlike plane- 
tary disk is (cf. Eq. 12 in lBaarsl[l973h : 



K 



1 — exp(- 



_ 1// 



with X = — 1 



(5) 



where / is called beam dilution factor. 

The apparent surface area of Uranus changed by 2.5% 
over the 20 years observing period discussed here. This is 
due to its oblateness, coupled with its high obliquity. We 
take this geometrical effect into account when discussing 
possible temporal changes of its disk averaged tempera- 
ture in Sec. 13. II The variations due to Neptune's oblate- 
ness are much smaller and ignored. 

3. Uranus and Neptune 

Because of their small apparent diameters of ~ 2", Uranus 
and Neptune are good calibration standards, at least for 
the larger and more sensitive mm-wavelength telescopes. 
In addition, they were never far apart in the sky during 



Table 1. Observations of the secondary calibrators at 
90 GHz for the years 1985 to 2005 and the adopted flux 
densities and sizes at this frequency. N is the number of 
observations. 



Source 


(Tl) 


rms 


iV 








[K] 


[%] 




[Jy] 


["1 




0.63 


8.2 


517 


3.77 


14 X 10 


NGC 7538 


0.41 


10.1 


351 


2.37 


20 X 18 


K3-50A(3) 


1.06 


8.7 


451 


6.00 


10 X 5 


NGC 7027'^' 


0.77 


7.9 


455 


4.58 


~ 10 



References: (1) ISteppe et all ljl993l l. for compariso n: (2) 
W3Q H: 4.08 ± 0.12 Jy, NGC 7027: 5.27 ± 0.15 Jy (|Ulicbl 

198111: f3) W30H:3 93 Jy, K3-50A: 6.31 Jy (Table 3 in 
Renter fc Krameil Il998l ); (4) a p proxim ate dimensions at 
mm- wavelengths f rom ISandelll (|l994l '). for comparison: 
iReuter fc Krameij l|l998h give at 90 GHz, W30H: 6" and 
K3-50A: 2". 
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Fig. 2. Ratio of Uranus and Neptune brightness temper- 
atures at 86-90 GHz normalized to the average ratio 1.075. 
The solid line is the result of a linear fit to the 81 observed 
data (small dots). Large dots are the corresponding yearly 
averages. 



the past 20 years (distance ^ 22°) and are therefore ideal 
for comparative me as ureme nts. 

Griffin fc Ort"CTl (Il993l) determined Uranus and 



Neptune brightness temperatures in the 0.35 to 3.3 mm 
region, using Mars as the primary standard. From mea- 
surements made in 1990 and 1992 they derive brightness 
temperatures of Uranus and Neptune at 90 GHz of 
rB(Uranus) = 136K and TB(Neptune) = 131 K. 

Since Uranus and Neptune are close together on the 
sky, a plot of the ratio of their brightness temperatures 
derived from simultaneous measurements made with the 
same telescope, may reveal also a relative change of their 
brightness temperatures. Such measurements of Uranus 
and Neptune made at frequencies between 86 and 90 GHz 
during 1986 and 2005, are shown in Fig. [2] These data 
were partly observed without complementary observations 
of the secondaries. However, only their relative ratios are 
discussed in this paragraph. 
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Fig. 3. Relative change of Uranus' brightness temperature 
at 90 GHz and at 8.6 GHz. The sohd hne is the resuh of 
a hnear fit to the 286 observed 90 GHz data (small dots). 
Their rms is 6%. Large solid dots are the yearly averages 
(cf. Table[2]). (Tb) is the average brightness temperature at 
90 GHz data of 134 K. Open circles show the 8.6 GHz data 
of KH06 normalized to Tut = 218 K which was extracted 
from Fig. 3 of KH06 for the period shown. 



Calculating the brightness temperature ratio from the 
ratio of antenna temperatures, we took into account the 
Uranus (U) and Neptune (N) apparent diameters at the 
time of the observation: rB(U)/TB(N) = TX(U)/TX(N) x 
/(U)//(N), with the beam dilution factor f{9s,eb) (cf. 
Eq.[5]). The average planetary diameter 0s was calculated 
from the equatorial and polar radii (JPL's HORIZONS sys- 
tem). In the case of Uranus, we took into account that 
the effective polar diameter changes with viewing geome- 
try (see below). 

The ratio is exp ected to be constant at 1.04 
( Griffin fc Ort(H]ll993[ ). The average ratio derived from 
the IRAM observations is 1.07, only slightly larger. We 
are, however, interested in the relative changes. The rms 
of the observed ratios is only 6%. 

A linear least-squares fit to the data 
( Bevington fc RobinsonI Il992l ) shown in Figure[2] re- 



sults in a slope of s = —0.4 10~^ per year, i.e. a drop 
by 8% between 1985 and 2005. However, the statistical 
error of the slope is very large and the linear-correlation 
coefficient is small (r = —0.37). In order to study where 
this variation may come from, we present in the following 
the Uranus and Neptune data independently. These data 
were observed near simultaneously with the secondary 
calibrators to derive their brightness temperatures. 



3.1. Uranus 

Uranus has a 84 year orbital period of which we have ob- 
served approximately one season (Fig. [3]), between solstice 
in the mid 1980s and almost equinox which will take 
place in 2007. Uranus is slightly oblate, the polar radius 
of 24973 k m is 2.4% smaller than its equatorial radius of 
25559km (|Lindallll987l . ll992f) . In addition, Uranus has an 
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Fig. 4. Relative change of Uranus' brightness temperature 
at 90 GHz and at 8.6 GHz as function of SEP latitude. 
Small filled dots show the individual 90 GHz measure- 
ments, large filled dots show the 90 GHz data binned into 
intervals of 15°. Errorbars are calculated from the scatter 
of the individual data. Open circles show the 8.6 GHz data 
of KH06. 



obliquity of 82°. In 1985, Uranus showed its south pole to 
the earth-bound observer, the sub-earth point (SEP) lati- 
tude was —82°, the apparent radius was nearly the equa- 
torial radius. It showed its equator in 2005 when the SEP 
latitude was —6° (cf. Fig.|4]). In summary, Uranus showed 
different regions and the disk area changed by ^ 2.5% over 
the past 20 years. The change of disk area has been taken 
into account in the following. 

The secular change in normalized brightness tempera- 
ture of Uranus between 1985 and 2005 at 90 GHz is shown 
in Fig. [21 These data show an overall scatter (rms) of, 
again, 6%. A linear least-squares fit to the data results 
in a slope of s = —0.48 10~^ per year, i.e. a drop by 9% 
between 1985 and 2005. Again, the statistical error of the 
slope is large and the linear-correlation coefficient is small 
(r = —0.35). The yearly averages vary between -f8% in 
1991 and -5% in 2005 (cf. also Table 3). Note that this 
drop is consistent within the errors with the change found 
when studying the ratio of Uranus and Neptune brightness 
temperatures (Fig. [5]), however with much fewer data. This 
indicates that Neptune temperatures are constant. 

Fi gurelSl also shows the results of iKlein fc Hofstadter 
(I2OO6) who studied Uranu s at 8.6 GHz (3.5 cm) (see also 
iHofstadter fc Butler] [2003h over the last 36 years and find 
significant variability. In late 1993, they find a strong 
temperature decrease, indicating a rapid, planetary-scale 
change. And indeed, the 90 GHz data also show a 
temperature-drop at this time. Here, the decrease is by 
-7%. 

A plot of the 90 GHz temperature variation against 
SEP latitude (Fig. [4]) shows an almost monotononic de- 
crease of almost 10%, from about -1-5% at SEP latitudes of 
about -85° to -4% at SEP latitudes of about -10°, indi- 
cating that Uranus brightness temperature depends on the 
viewing angle. This decrease is detected at the ^ 2a-level 
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Fig. 5. Uranus thermal profile P{T) (jLindall [l987h and 

vertical distribution of NH3 and CH4 mixing ratios. The 
values of the NH3 and CH4 mixing ratios below their 
condensation levels are equal to 6xl0~® and 2.3x10"^, 
respectively. These have been used as input to a radia- 
tive transfer model to calculate the weighting functions of 
Uranus' con tinuum. The model had previously been used 
for Jupiter (jMoreno et al.l 120011 ). The weighting function 
at 8.5, 90 and 230 GHz are shown in black. 



Table 2. Yearly averages and rms of Uranus' brightness 
temperature (cf. Fig[3]). Tq is the yearly average, (Tb) is 
the overall average of 134 K. is the number of observa- 
tions. 



Year 


Tb/{Tb) 


rms 

[%] 


iV 


1986 


1.05 


4.6 


8 


1987 


1.04 


6.1 


5 


1991 


1.08 


2.7 


6 


1992 


1.06 


2.6 


12 


1993 


1.06 


6.7 


14 


1994 


0.99 


3.9 


10 


1995 


0.98 


7.1 


11 


1996 


0.99 


6.6 


14 


1997 


0.99 


7.3 


16 


1998 


1.00 


6.4 


18 


1999 


1.03 


5.0 


19 


2000 


1.03 


7.4 


21 


2001 


1.00 


3.7 


38 


2002 


0.97 


4.0 


32 


2003 


0.97 


4.5 


19 


2004 


0.97 


4.7 


28 


2005 


0.95 


7.2 


14 



as the binned data have an average rms of 5.4% (Table[3l). 
On average, the 8.6 GHz data of KH06 show a similar de- 
pendence with year and SEP latitude, however with larger 
scatter. 

FigureO shows that higher frequencies probe succes- 
sively higher layers of the Uranus atmosphere. The con- 
tinuum at 8.5 GHz stems from the high pressure zone, i.e. 
11.0 to 5.7 bar at low altitudes, where the opacity is sen- 
sitive to the NH3 vertical distribution and the thermal 
structure. In contrast, 90 GHz emission probes the region 
of pressures between 6.1 and 4.7 bar, higher in the at- 
mosphere, where the opacity is also a function of NH3, 
but much more sensitive to the thermal structure. At still 
higher frequencies, the 230 GHz continuum emission arises 
from still higher atmospheric levels (3.5-0.8 bar), near the 
tropopause. The continuum at 230 GHz is mainly sensitive 
to the thermal structure, and only slightly to the CH4 ver- 
tical distribution, but not to NH3. 

It is important to note that between 8.5 and 90 GHz, 
the difference of the opacity sensitivity to NH3 and to the 
thermal structure, may explain the differences of the rel- 
ative changes of Uranus' brightness temperature between 
8.5 and 90 GHz shown in Figures[3] and [H 

3.2. Neptune 

Figureini shows the brightness temperatures of Neptune, 
normalized to the average temperature of (Tb) of 128.3 K. 

^ Using the equatorial radius for the calculation of the disk 
area at all times, i.e. not taking into account the change of disk 
area with viewing angle, would lead to a slightly steeper slope 
and a decrease of —7% at SEP latitudes of —10°. 



Table 3. Normalized Uranus brightness temperatures as 
function of SEP latitudes binned into intervals of 15° 
width (see also Fig.H]). N is the number of observations. 



SEP 


(Year) 


Tb/(Tb) 


rms 


N 


Latitude 




[K] 


[%] 




-85.0 


1986.70 


1.05 


5.0 


13 


-70.0 


1991.92 


1.06 


3.6 


13 


-55.0 


1994.16 


1.03 


6.7 


44 


-40.0 


1998.17 


1.00 


6.5 


57 


-25.0 


2001.67 


0.99 


5.5 


108 


-10.0 


2004.55 


0.96 


5.6 


51 



A linear least-squares fit indicates a very small drop of 
only 2% between 1985 and 2005. In addition, the statisti- 
cal uncertainty of the fitted slope is large, 80%. We con- 
clude that Neptune temperatures are constant within the 
8% observational error. 

4. Summary 

By careful selection of data from the archived pointing 
measurements made at the IRAM 30m telescope during 
the period 1985 to 2005, we have been able to study the 
long-term behaviour of Uranus and Neptune at 90 GHz. 
The large number of selected observations allowed a sta- 
tistically meaningful analysis. 

1. We obtained 286 observations of Uranus. The scatter 
of the derived normalized brightness temperatures is 
6%. 
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Year 

Fig. 6. Normalized Neptune brightness temperature. 
Small dots and crrorbars show the individual observations 
and their rms errors. Dashed lines show the 8% (10 K) rms 
scatter of all 86 observations. 

We find a systematic variation of the Uranus bright- 
ness temperature of ~ 10% with sub-earth point (SEP) 
latitude when the orientation of Uranus changes from 
south-pole view to equator view. This effect is detected 
at the ~ 2(7 level, the mean scatter of the data averaged 
in 15° bins is 5%. A similar change is indicated by the 
8.6 GHz observations of KH06. The year 1993 shows a 
rapid decrease of brightness temperatures at 90 GHz of 
~ 7%, which corresponds to a similar decrease seen by 
KH06, indicating a temporal change of a large fraction 
of the Uranus atmosphere. More observations at other 
frequencies and improved models of the Uranus atmo- 
sphere, taking into account its latitude structure, are 
needed to better determine the variations with SEP. If 
Uranus is to be used as calibration source at 90 GHz, 
the latitude dependence of its brightness temperature 
needs to be taken into account. 
2. For Neptune, we only obtained 71 observations tied 
to secondary calibrators and 81 simultaneous obser- 
vations with Uranus. Both show that its brightness 
temperature stays constant at a level of better than 
8%. 
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